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ABSTRACT 

AA Tau, a classical T Tauri star in the Taurus cloud, has been the subject of intensive photometric 
monitoring for more than two decades due to its quasi-cyclic variation in optical brightness. Beginning 
in 2011, AA Tau showed another peculiar variation - its median optical though near-lR flux dimmed 
significantly, a drop consistent with a 4-mag increase in visual extinction. It has stayed in the faint 
state since. Here we present 4.7 /rm CO rovibrational spectra of AA Tau over eight epochs, covering 
an eleven-year time span, that reveal enhanced ^^CO and ^^CO absorption features in the Jiow ^ 13 
transitions after the dimming. These newly appeared absorptions require molecular gas along the line 
of sight with T~500K and a column density of log (Ni 2 co) ~18.5cm“^, with line centers that show a 
constant 6 km s“^ redshift. The properties of the molecular gas confirm an origin in the circumstellar 
material. We suggest that the dimming and absorption are caused by gas and dust lifted to large 
heights by a magnetic buoyancy instability. This material is now propagating inward, and on reaching 
the star within a few years will be observed as an accretion outburst. 

Subject headings: protoplanetary disks -stars: individual (AA Tau)-stars: pre-main sequence - stars: 
variables: T Tauri, Herbig Ae/Be 


1. INTRODUCTION 

Sun-like stars are born surrounded by disks which reg¬ 
ulate the angular momentum redistribution required for 
accretion. It is also in these disks that planets form. 
Understanding the angular momentum transport and 
accretion processes in circumstellar disks is thus crit¬ 
ical for our understand ing of star and planet forma¬ 
tion (| Turner et al.ll20f^ . Photometric and spectroscopic 
monitoring of young stars and their accretion disks is an 
important way to study the dynamics at play. 

The classical T Tauri star AA Tau has provided 
many clues to the magnetospheric accretion process 
that channels material from the inne r edge of the disk 
onto the central star (|Bouvier et al.l 1199911 . A A Tau 
shows a quasi-cyclic variation with a period of 8.5 days 
in its optical photometry with an ampl i tude of 1.4 
mag i n the BVRI band s (iBouvier et all 119991 120031 
l2007t iMenard et al.l l2003t iGrankin et al.l 1200711 . This 
peculiar behavior is interpreted as a result of peri¬ 
odic occultations of the star by a non-axisymmetric in¬ 
ner disk warp, perhaps driven by a misalignment be- 
tween the stellar magnetic dipole and rotation axes 
(IBouvier et al.lflO^ l2003HTerauem fc PaDaloizoull2000t 

O’Sullivan et al.l 120051 : iDonati et al.l 120101 : lEsau et al.l 
201411 . Such occultations require that the AA Tau disk 
has a nearly edge-on (z ^ 71 — 75°) orientation, which 
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has indeed been seen in optical scattered light images 
(|Cox et al.ll20I3h . 

Recently, another peculiar variation was discovered 
(|Bouvier et al.l [20I3I1 . In 2011, AA Tau dimmed ^ 2 
- 4 mag in the V band after maintaining a U-band 
magnitude between 12.3 and 14.3 mag for m ore than 
two decades (1987-2010, iCrankin et al.l 120071 see Fig¬ 
ure [IJ . Its visual extinction increased by 4 mag based 
on the JHK color changes, and A A Tau has remained 
in the faint state since (September 20 14, current V- 
band mag=14-17, IBouvier et al.l I2013L and Grankin, 
priv. commun.). Surprisingly, no significant change was 
found in the mass accretion rate onto the central star 
by comparing accretion tracers observed before and one 
year after the d imming started (ba sed on Ha, H/3 and 
[Hel] line fluxes. IBouvier et al.ll2013l) . This is in contrast 
to other observations of disk variability, whi ch are usu¬ 
ally a cc ompanied by a c cretio n rate changes (|Sitko et al.l 
l2012fl . IBouvier et al.l (l2013fl concluded that the sud¬ 
den dimming is due to additional circumstellar extinc¬ 
tion along the line of sight and suggested that the extra 
extinction is possibly produced by a non-axisymmetric 
overdense region in the outer disk which recently moved 
into the line of sight via disk Keplerian rotation. Despite 
the importance of discovering the origin of the enhanced 
extinction, photometric data alone provide limited in¬ 
sight into the physical and dynamical properties of the 
extra material along the line of sight, so we turn to spec¬ 
troscopy to measure gas kinematics and temperature. 

Here we present new observational constraints on the 
origin of the dimming of AA Tau in the form of eight 
epochs of CO 4.7 /rm spectra taken over an eleven-year 
time span. Significantly enhanced ^^CO and ^^CO ab¬ 
sorption components are seen in the spectra taken after 
the dimming started, at a velocity redshifted with respect 
to the star. The spectra can be used to measure the prop- 
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erties of the absorbing gas, providing new constraints to 
the origin of both the dimming and absorption. Indeed, 
since the optical extinction and CO absorption appear 
at similar times, and have persisted since, the absorbing 
gas and obscuring dust most likely arise from the same 
event (s). We argue that a non-axisymmetric region in 
the Keplerian disk is unlikely to account for the extra 
absorption and velocity shift, and suggest that the en¬ 
hanced extinction and extra absorption may be produced 
by a disk instability-driven infalling flow. 

2. OBSERVATIONS 

We carried out high resolution spectroscopic obser¬ 
vations of AA Tan at eight different epochs between 
November 2003 and September 2014. Figure 1 shows 
the dates of the M-band observations, overlaid on the 
V-band photometry time series of AA Tau over the 
last decade. A detailed observation log is provided in 
Table K_ Most of our spec tra were taken with NIR- 
SPEC ( McLean et al.l 119981 1. a high resolution spec¬ 
trometer (R^25,000) at the Keck II telescope, as part 
of a large NIRSPEC survey of protopl anetary disks 
(|Blake fc BoogertI 12004 iSalvk et ‘ahl 1200911 . in both the 
native seeing and adaptive optics (AO) modes. The 
M-band (4.7 /im) echelle spectra were observed with a 
0."43x24"slit (native seeing observations) or a 0."041 x 
2."26 slit (with AO). Two spectral settings were used, 
covering wavelengths between 4.65 and 5.15 /rm, with 
gaps between orders. The total wavelength range covers 
a large portion of the v = (1-0) fundamental rovibrational 
band of ^^CO, i.e., the first two R-branch lines (the total 
angular momentum quantum number of the lower energy 
state Jiow=0,l) and the low/mid P-branch (Jiow=l-12 
and Jiow=30-40), as well as an accretion tracer, the Hi 
Pf/3 transition. The AA Tau M-band spectrum of Octo¬ 
ber 2007 is obtained from archival VLT CRIRES data, 
taken as part of a large survey on protos t ars and pro¬ 
topla netary disks (|PontoDDidan et al1l20lH iBrown et al.l 
1201311 . The CRIRES spectrum covers a wavelength range 
between 4.65 and 4.77 fim at a spectral resolution of 
-95,000. 

Objects were observed in nod pairs, with subsequent 
pairs subtracted from one another to remove telluric 
emission features. The differenced images were flat field 
corrected and then averaged to increase the signal-to- 
noise ratio. We extracted 1-D spectra from t he 2-D 
avera ged images using the optimal extraction (|Hornel 
1198611 . The wavelength calibration was derived by fit¬ 
ting a fourth order polynomial to the telluric emission 
lines within the same echelle order. Telluric absorption 
features were removed by dividing the AA Tau spectra by 
that of HR 1620 (A7 spectral type), taken close in time 
and airmass. To investigate the continuum flux varia¬ 
tion with time, we compared the raw photon counts of 
AA Tau with those of HR 1620 in the continuum around 
4.7 ^m at different epochs. We found the flux ratio varied 
less than 30% from the midpoint of all values obtained 
over the past eleven years, with no apparent correlation 
between the continuum flux and CO emission line inten¬ 
sities. We thus assumed the 4.7 ^m continuum is con¬ 
stant and flux calibrated the 1-D spectra by normalizing 
the continuum to the interpolated IRAC fluxes between 
3.6 Ai m (352 mJy) and 4.5 /rm (332 mJvl dLuhman et al.l 
l2006f) . It is worth noting that we do not expect to detect 


Table 1 

AA Tau Observation Log 


Obs. Date 
(UT) 

Spectral Range 
{^lm) 

Int. time 
(min) 

Instrument 

AO“ 

2003-11-03 

4.65-4.72, 4.95-5.03 

12 

NIRSPEC 

N 

2004-12-27 

4.65-4.72, 4.95-5.03 

12 

NIRSPEC 

N 

2004-12-30 

4.70-4.78, 5.02-5.09 

8 

NIRSPEC 

N 

2007-10-12 

4.65-4.77 

24*’ 

CRIRES 

Y 

2010-12-15 

4.70-4.78, 5.02-5.09 

64 

NIRSPEC 

Y 

2012-12-31 

4.65-4.72, 4.95-5.03 

32 

NIRSPEC 

Y 

2013-01-01 

4.70-4.78, 5.02-5.09 

24 

NIRSPEC 

Y 

2013-10-19 

4.65-4.72, 4.95-5.03 

24 

NIRSPEC 

Y 

2013-12-23 

4.65-4.72, 4.95-5.03 

24 

NIRSPEC 

Y 

2014-09-06 

4.65-4.72, 4.95-5.03 

47 

NIRSPEC 

N 


Note. — a: Observations with adaptive optics; b: The final spec¬ 
trum combines observations with two partially overlapping spectral 
settings with 10 and 14 min. of integration time, respectively. 


continuum variability caused by the enhanced extinction, 
because an Ay of 2-4 mag produces a flux change of 
less than 10% at 4.7 /im using a Mathis extinction law 
(IMathislll99nll . 
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Figure 1. V-ban d photometry time series of AA Tau 
JBouvier et al.l [20131) overlapped with observational dates of CO 
M band spectra (red vertical dash lines). The 2013 December 
photometry data are from Grankin and Bouvier (priv. commun.). 


3. RESULTS 

Figured] displays our CO M-band spectra of A A Tau 
taken between November 2003 and September 2014 (data 
obtained within three days were consistent and merged 
into a single spectrum). Readily detectable v = (1-0) fun¬ 
damental rovibrational ^^CO emission lines are seen in all 
eight epochs. The ^^CO emission components appear to 
be double-peaked and spectrally resolved with an average 
FWHM of —130km s“^. The double-peaked line shape 
is expected for emission from a highly inclined, Keplerian 
disk. By fitting the CO line profile, ISalvk etldl (|2011[1 
suggested the molecular emission arises from —O.IAU. 
Our AO data show no CO spectro-astrometric signatures 
along the disk major axis to a 3 (t limit of 0.19 AU, consis¬ 
tent with the previous result. In five of the eight epochs, 
the lines were asymmetric, with the blue-shifted peak 
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stronger than the red-shifted peak (Figure IS])- The line- 
to-continuum ratio also varied significantly, from 0.15 to 
0.38, but there is no apparent link between the emission 
line variability and the F-band extinction change. 

The most striking discovery is the emergence of deep, 
low- to moderate-J (Jiow<13) and absorp¬ 

tion below the continuum in spectra taken after the 
optical dimming (see Figure [2] and [3]). The CO ab¬ 
sorption components have a FWHM^14.2km s“^, a 
value close to the width of the nominal instrumental 
profile (FWHM^12km s“^). The mid-J CO lines be¬ 
tween 4.95-5.09 /rm (of lower state excitation energies 
Fliow ^2000 K), however, did not show enhanced absorp¬ 
tion after the dimming, suggesting the absorbing gas can¬ 
not be hotter than ^1000 K. For all the spectra taken be¬ 
fore the dimming, the center of the double-peaked emis¬ 
sion line profile aligned well with the stellar velocity (he¬ 
liocentric velocity V* = 16.1±2km s~^. lHartmann et al.l 
1198611 . In contrast, the absorption line centers from 2011 
and after show a clear and constant redshift of ^6±1 km 
s“^. 

Because the goal of this paper is to shed light on the 
origin of the dimming and absorption towards AA Tau, 
we focus here on the properties of CO absorption lines. 
The varying CO emission of A A Tau will be discussed 
in a separate paper which investigates the variation and 
asymmetry of M-band CO emission in a large sample of 
disk sources (Crockett et al. 2015, in prep). We have 
taken additional L-band NIRSPEC observations before 
and after the dimming to search for additional molecules 
such as water in the disk gas; but the complex pho- 
tospheric structure present at these wavelengths means 
that only very deep absorption lines would be recognized, 
and no such features have been found. 

3.1. The physical properties of the absorbing gas 

The line profiles of the absorption components are pos¬ 
sibly varying (see Figure|3]), but variations in the Earth’s 
atmosphere can also mimic the line profile change since 
both the telluric and AA Tau features are marginally 
resolved, at best. Nevertheless, the line fluxes are less 
affected and the rotational ladder of CO v = (1-0) lines 
provides abundant information on the physical proper¬ 
ties of the absorbing gas. Under the assumption of a 
Boltzmann distribution for the rotational states, the line 
flux ratios provide a direct constraint to the excitation 
temperature of the gas. Once the excitation tempera¬ 
ture is known, the total column density of the absorbing 
gas can be derived from the absolute strengths of the 
absorption components. 

This assumption of local thermodynamic equilibrium 
is reasonable in this case because the CO?; = 0,J^13 
states are expected to be thermalized due to their mod¬ 
erate critical densities of ^ 10^ cm“^. Given the lack 
of constraints on the geometry of the absorbing gas, we 
make the simplest assumptions, i.e. the gas is isothermal 
and has a constant number density along the line of sight. 
Our model thus has three free parameters: the gas tem¬ 
perature T, the column density N along the line of sight, 
and the intrinsic line width a (needed because the ^^CO 
and ^^CO absorption lines are largely spectrally unre¬ 
solved at R=25,000). The line width is an upper limit to 
the local turbulence velocity since the absorbing gas may 


exist over a range of radial distances and azimuthal an¬ 
gles. The differential projected radial velocity from gas 
at different locations can broaden the line profile. 

We measured the equivalent width {W) of each ab¬ 
sorption line asW = J (/cont — d(u))//cont dv. The ^^CO 
equivalent line widths provide a particular challenge since 
the absorption is built on top of emission lines. Because 
our goal is the study of the additional absorption that 
appears after the dimming, we measured the equivalent 
widths of the ^^CO absorption from the difference be¬ 
tween spectra, specifically those from December 2004 and 
January 2013 because they cover the largest wavelength 
range and thus provide the best constraints on the exci¬ 
tation temperature. Each absorption component in the 
residual is then fit with a gaussian function rather than 
directly integrating over the line profile since certain ve¬ 
locities are missing in regions that overlap with saturated 
atmospheric lines. The gaussian fits are then divided 
by the continuum flux to calculate the ^^CO equivalent 
widths, while those for ^^CO were measured directly from 
continuum subtracted January 2013 spectra. 

The theoretical equivalent width of a given transition 
is calculated using ly = Icontexp{—sNj(j)y), where Nj is 
the CO column density of a rotational state J, (f>y is the 
intrinsic line profile, and s is the integrated cross section 
(cm^ s“^), or 


s = 


gu . 

87ri/2 gi 



( 1 ) 


The spontaneous decay rate, A^i is taken from HITRAN 
(Rothman 2005), while a Voigt function is used for (fy. 
The Doppler width of Voigt function, = \/2cj^ con¬ 
tains a as the third free parameter. We used an average 
ISM value of ^^CO/^^CO number ratio of 69 to scale 
the c olumn densities of ^^CO to those of ^^CO (IWilsonI 

[IMl)- 

The fit results are summarized in Eigure01 Our best fit 
has T = 494lgg^ K, with a total column density of ^^CO 
along the light of sight of log(Ni 2 co) = 18.5±0.1cm“^ 
and an intrinsic width a of 2.2±0.1 kms“^. We stress 
that this result is not very sensitive to the input iso¬ 
topic ratio; changing the ^^CO/ ^^CO value from 69 to 
30 would result in T = 427K, log(Ni 2 (;;;o) = 18.1, and 
cr = 2.6. The best-fit column density decreases in this 
case because the ^^CO absorption lines are saturated and 
only the optically thin ^^CO lines provide the main con¬ 
straints on the column density. A smaller ^^CO/ ^^CO 
ratio thus produces a smaller N(^^CO). 

For the remaining epochs, we did not fit the three 
parameters (T, N and cr) simultaneously because the 
available data often cover a smaller wavelength range, 
with many absorption lines that are heavily contami¬ 
nated by saturated atmospheric absorption due to unfa¬ 
vorable Earth-induced Doppler shifts. Since the gas col¬ 
umn density is most sensitive to the ^^CO line strengths, 
we derived the CO gas column density in these epochs 
by fixing the T and cr parameters with the best-fit values 
from January 2013. In the three 2013 spectra the derived 
column densities of the absorbing gas vary by <30%. Eor 
the September 2014 epoch, the column density decreased 
to log (N^^CO) = 18.2±0.1 cm“^, about a half of that 
in 2013. In the pre-dimming epochs, the October 2007 
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Figure 2. CO M-band spectra of AA Tau before and after the 2011 F-band dimming, (a) CO rovibrational spectra between 4.64 and 
4.78 fim. The blank regions represent regions with telluric transmission <25%, which we have masked out due to high noise levels. The 
rest wavelengths of the ^^CO transitions are marked as small vertical lines at the top of the panel. The region highlighted by the dashed 
box is enlarged in panel (b). (b) A close up of the 4.64-4.73 /rm region, highlighting the ^^CO absorption lines that appear after 2011. 
Rest wavelengths of the ^^^CO v= 1-0 lines are marked with dashed lines, (c) High-J CO spectra between 4.95 and 5.09 /rm. 
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Figure 3. Time variation in the low-J and line 

shape(s). The average line shape of each epoch combines the 

V =(1,0) Jiow ^12 lines, that for averages the isolated 

from R12 and 13 lines. There is no line profile for 

December 2010 epoch because the spectra do not include the two 
transitions. All of the lines are normalized to the continuum 
fiux. The two vertical dashed lines indicate the locations of stellar 
velocity and a red shift of 6 km s“^. The horizontal dotted lines 
denote the onset of dimming in 2011. 



log [cm ^] = 18.5±0.1 



Figure 4. Top: Best model of the equivalent widths of the 
and absorption components in January 2013. Error bars 

indicate the 1 cr uncertainties, filled circles (black ones are 

data points and red circles are best model fitting results), filled 
squares (the data points are in black while the model results 

are in blue). Bottom: Probability distribution surface contours of 
T and Ni 2 qo, at 68%, 95%, 99% confidence levels. Marginalized 
one dimensional posteriors are shown on each side. The shaded 
regions depict the 68% confidence areas. 


spectrum provides the most stringent column density up- 
per limits, suggesting log (N^^CO) < 17.5 cm“^ at 99% 
confidence (T and a are fixed). Because the lines 

are optically thin, these column density changes can be 
directly seen in Figure [3l 

The contemporary existence of the absorbing gas and 
the obscuring dust suggests that they are likely to arise 
from the same event (s). Assuming the gas and dust are 
co-spatial, we can derive an average gas-to-dust mass ra¬ 
tio in the absorbing material along the line of sight, by 
combining the column density with the increase in 

visual extinction associated with the dimming. Assum¬ 
ing a maximum fractional abundance of nco/?T-H=10~^ 
in protoplanetary disks, the total column density of ab¬ 
sorbing gas is at least 3.2x10^^ cm“^, resulting in a 
Vn/Ay >8x10^^ cm“^ mag“^ - about a factor of 4 
larger than the average ratio measured in nearby molec¬ 
ular clouds (Vuong et al 2003), suggesting a gas-to-dust 
ratio of ^400 in the absorbing gas. Such an increased 
gas-to-dust mass ratio is expected if the absorbing gas at 
large disk heights has experienced significant dust grain 
growth and/or settling towards the mid-plane. The slope 


of the spectral energy distribution of AA Tau in (sub)mm 
wavelength range is co nsistent with the expectatio n of 
significant dust growth ([Andrews fc Williamsll2007ll . 

4. THE ORIGIN OF THE ABSORBING GAS 

The fairly sudden emergence of enhanced extinction 
and absorption indicates that some obscuring material 
recently appeared or moved into our line of sight. We ar¬ 
gue that this new material is unlikely to be a foreground 
source. Eirst, the absorbing gas is too warm 500 K) 
compared to the typical interst ellar diffuse cloud exc i- 
tation t emperatures of 10-100 K (|Snow fc McGallir200fill . 
Second, iBouvier et al.l (1201311 found that the enhanced 
extinction was restricted to AA Tau since three field stars 
nearby (projected distance between 26"and 53") did not 
exhibit any luminosity variation. This projected distance 
is much smaller than that of a typical gas clump in molec¬ 
ular clouds at 140 pc. Eor example, the Jeans length of 
gas with T of 500 K and Nhj of 10"^ cm“^ is 1.3 pc, i.e., 
a projected distance of 1915". 

The absorbing gas is also unlikely to be in circumstellar 
envelope surrounding AA Tau because various evolution- 
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ary indicators suggest that AA Tau is a Class II object, 
i.e., a disk with little or no envelope (|Robitaille et alJ 
l2007f) . Moreover, gas in a c ircumstellar envelope shou ld 
be much cooler than 500 K (jvan Dishoeck et al.ll201,‘l|l . 

The final possibility is that the absorbing material lies 
in the disk. We can place some constraints on its loca¬ 
tion based on the gas temperature of ^500 K. Utilizing 
existing thermo-chemical models in the literature which 
are constructed for classical T Tauri disks similar to AA 
Tau, we can use this temperature to estimate the vertical 
and radial location of the absorbing gas. In the vertical 
direction, due to the high inclination of the disk and the 
optical depth of the infrared CO transitions observed, 
the absorbing gas must be high in the disk atmosphere 
but deep enough that sufficient CO column density is 
present in order to shield itsel f from photodissoc iation 
by high-energy stellar photons (jVisser et al.ll200^ . The 
dominant heating sources for this warm molecular layer 
are X-ray and UV radiation (e.g. iGlasseold et al.ll2004 
iGorti et al.ll20HTl . The X-ray luminosity of AA Tau is 
~ 10 ^° erg s“^ in quiescence and ^ 10 ^^ erg s~^ dur¬ 
ing b rief flares (|Schmitt fc Robra'd^ l2007t iGrosso et al.l 
l2007fl , all wit hin the typical X-ray l u minosity of classica l 
T Tauri stars (|Feigelson et al.ll^OOSH . iNaiita et al.l (120111) 
computed a thermo-chemical model for classical T Tauri 
stars with an Lx of 10^° erg s“^. Their model shows that 
CO starts to survive beyond a vertical column density of 
Nh 10^^cm“^. At such columns, a temperature of 
500 K requires a radi us that is less than 10 AU distant 
from the central star. iWalsh et al.l (|2010l I 2012 H . for ex¬ 
ample, used the X-ray, UV continuum and Lya spectra 
of TW Hya in their thermal chemical model, and their 
calculations predict that the high altitude CO molecular 
layer inside ^10 AU is warmer than 500 K. In summary, 
the warm temperature of the absorbing gas suggests a 
location of ^10 AU in radius. 

We conclude with a discussion of the possible origins 
of the suddenly enhanced extinction and CO absorptions 
in the AA Tau disk: 

( 1 ) A non-axisymmetric overdense region in the disk 

iBouvier et al.l (1201311 suggested that the enhanced ex¬ 
tinction is caused by a non-axisymmetric structure that 
recently moved into our line of sight due to Keplerian ro¬ 
tation. However, we argue that this scenario becomes 
less likely considering the CO absorption components 
presented here. Because the M-band continuum and 
the CO emission both arise from very small regions 
(R<0.5AU), only gas and dust located inside a small 
solid angle can produce the enhanced extinction and ab¬ 
sorption. Therefore, the absorption lines should be cen¬ 
tered at the stellar velocity. However, the CO absorption 
lines display a constant redshift of ^6 km s“^ over the 
past two years. This redshift cannot be explained by the 
uncertainty of the stellar velocity of AA Tau as all pre¬ 
vious measurem ents of the stellar veloc ity are consistent 
within 2km s ~^ ([Hartmann et al.lIlQS^ of 16.1±2km s“^, 
Bouyier et al.ll200^ of 17.1±0.9km s“^ and iDonati et all 
20ir)l of 17.2T0.1 km s“^). A highly eccentric orbit might 
be able to explain the redshift, but the CO emission lines 
are well centered at the stellar velocity - suggesting at 
least the disk gas in the CO emitting region is in circular 
motion. If the non-axisymmetric structure rotates at Ke¬ 
plerian speeds, its orbital period is then at least 27 years. 


the total time span of AA Tau in the bright state and 
the duration of the dimming. Thus the structure must 
lie beyond >8.4 AU (assuming M*=O. 8 M 0 , Bouvier et 
al. 1999). 

(2) Inward flow driven by disk instability 

Another possibility is the absorbing material has been 
newly lifted to large heights due to an instability of the 
disk’s magnetic fields, leading to an accretion outburst 
involving the inward transport of gas and dust. The 
buoyant rise of the magnetic fields generated in magneto- 
rotational turbulence has been demonstrated in stratified 
isothe rmal shearing-box m agnetohydrodynmic calcula¬ 
tions (iMiller fc Stond [2000[) and in calculations includ¬ 
ing the stabilizing effect of the external starlight heat¬ 
ing (iHirose fc Turned 1201111 . In all cases, the magnetic 
fields dominate the pressure above about two density 
scale heights. The fields are erratic, with stronger activ¬ 
ity coming in bursts, if the magnetic diffusivity is near 
the thres hold for switching off the magneto-r otational 
dynamo (|Turner et al.ll200A ISimon et al. I2011D . 

As such material propagat es inward, an acc retion out¬ 
burst can potentially result (|Zhu et al.l 1201011 . Episodic 
accretion outbursts have been observed in some pre-main 
sequence low-mass stars (e.g. see a review by Audard 
et al. 2014). Well-known examples include FU Orionis 
objects (with typical accretion rates during an outburst 
of 10 “’^ Mq yr“^ and a decay timescale of ~100 years) 
and EX Lupi objects (smaller but more frequent out¬ 
bursts with a typical outburst accretion rate of 10 “^ Mq 
yr“^). Gould AA Tau be a pre-FUor/EXor? AA Tau 
has a quiescent accretion rate between 10 “® and 10 “® 
Mq yr“^, and has shown n o significant accreti on rate 
increase since the dimming (IBouvier et al.ll^OlSD . Sim¬ 
ilarly, we did not detect significant changes in the ac¬ 
cretion tracer Hi Pf/3 around 4.65 /rm before and after 
the absorption appeared. In Section 3.1, we measured a 
column density in the extra gas component of 3.2x10^^ 
cm“^. Assuming this gas is present in a wall at 10 AU 
with a scale height of 1 AU, and that this material will fall 
onto the star within three years, the accretion rate would 
rise to ^1.3x10“^ Mq yr“^. Based on the evolutionary 
stage of AA Tau and our estimate of the accretion rate 
of the current outburst, we suggest A A Tau is unlikely to 
be a pre-FUor object but could be a pre-EXor candidate. 

There is one additional hint that suggests episodic 
outbursts have p reviously occurred in AA Tau, in that 
ICox et all (I2013f) found a chain of Herbig-Haro knots 
above the AA Tau disk in [SH] images. Such Herbig-Haro 
objects are usually a ssociated with elevated a ccretion ac¬ 
tivity in young stars (iReinurth fc Ball's^l200lD . The spac¬ 
ing of Herbig-Haro knots in the A A Tau [SH] images, 
together with the inferred speed of knots >200 km s“^, 
suggests a time between ejections of at most a few years, 
which appears consistent with the smaller and more fre¬ 
quent outbursts in the EX Lupi type. 

The observed redshift in the absorption suggests the 
molecular gas is propagating inward. The magnitude of 
the speed is, however, puzzling: 6 km s“^ is too high 
for a viscous disk but too low for a free-fall velocity. 
With a central star of 0.8M 0 , the free-fall velocity 
^J2G^LflJr is 38 km s“^ at 1 AU. A disk with viscosity 
V will evolve at radius r on a timescale of t^ ^ r^/iz. 
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We adopt the usual parameterization of the viscosity 
as ^ = aCs/fi (cs is the sound speed, the Keplerian 
angular v elocity) and use the a ~ 0.1 found in FU Ori 
outburst (|Zhu et al.l 1200711 to estimate a viscous disk 
timescale. The radial velocity so obtained is rjt^ 


akT 

firuH 

s-i. 


v'r/GM*=0.07x 


0.1 500 K\AU J 


( r A/2, M, n-1/2 

\AU) Io.SMqI 


The instability driven inward flow scenario is spec¬ 
ulative since certain critical disk physical parameters, 
such as the magnetic field and viscosity, are largely un¬ 
known. Monitoring the physical properties of the absorp¬ 
tion components and the accretion rate over the next few 
years will provide crucial constraints on the origin of the 
optical dimming. With a radial velocity of 6 km s“^, mat¬ 
ter within 10 AU from the central star (T~500K) should 
drive an accretion rate increase over the next few years if 
the absorption lines are caused by infalling gas. Any sub¬ 
sequent outburst’s accretion rate and duration combined 
would yield the amount of material involved, while the 
duration of the absorption would suggest a location for 
the triggering instability. Furthermore, assuming the in¬ 
flowing material rotates at near-Keplerian velocities, the 
duration of the absorption constrains the azimuthal ex¬ 
tent of the absorber. For example, at a distance of 5 AU, 
the three-year ongoing absorption suggests an angular 
extent of at least 90 degrees. In short, the dimming and 
absorption observed in AA Tau may contain important 
clues into the transport processes in disks. 
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